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ABSTRACT 

Information about stellar magnetic field topologies is obtained primarily from high-resolution circular polarisation (Stokes V) obser¬ 
vations. Due to their generally complex morphologies, the stellar Stokes V profiles are usually interpreted with elaborate inversion 
techniques such as Zeeman Doppler imaging (ZDI). Here we further develop a new method of interpretation of circular polarisation 
signatures in spectral lines using cumulative Stokes V profiles (anti-derivative of Stokes V). This method is complimentary to ZDI and 
can be applied for validation of the inversion results or when the available observational data are insufficient for an inversion. Based 
on the rigorous treatment of polarised line formation in the weak-field regime, we show that, for rapidly rotating stars, the cumulative 
Stokes V profiles contain information about the spatially resolved longitudinal magnetic field density. Rotational modulation of these 
profiles can be employed for a simple, qualitative characterisation of the stellar magnetic field topologies. We apply this diagnostic 
method to the archival observations of the weak-line T Tauri star V410 Tau and Bp He-strong star HD 37776. We show that the 
magnetic field in V410 Tau is dominated by an azimuthal component, in agreement with the ZDI map that we recover from the same 
data set. For HD 37776 the cumulative Stokes V profile variation indicates the presence of multiple regions of positive and negative 
field polarity. This behaviour agrees with the ZDI results but contradicts the popular hypothesis that the magnetic field of this star is 
dominated by an axisymmetric quadrupolar component. 
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1. Introduction 

Stellar magnetism represents an important, though not fully un¬ 
derstood and poorly constrained, ingredient of the theories of 
stellar formation and evolution. Although magnetic fields are 
believed to play an important role in many astrophysical situ¬ 
ations, their direct detection and characterisation is often very 
challenging. The stellar surface magnetic fields with typical 
strengths ranging from a few Gauss to several kilo-Gauss are 
commonly diagnosed with the help of the Zeeman effect. The 
broadening and splitting of magnetically sensitive spectral lines 
lines becomes apparent in high-resolution spectra when the field 
strength exceeds ~ 1 kG, allowing one to measure the mean field 
modulus ([Mathys et al.| 1997 [ [Kochukhov et al.|2()0^ [Reiners & 


Basri 


& Va] 


2Q07|) and identify different field components (| Johns-Krull 

enti|1996t |Shulyak et al.|2014] ) in strongly magnetised ob 


lar polarisation (Stokes V) signatures. Several integral magnetic 
observables, fo r example the mean longitudin al magnetic field 
( [Mathys I [l 99 1| ) and crossover ( [Mathys [1995] ), can be derived 
by computing the wavelength moments of such simple Stokes 
V profiles. These observables are interpreted by fitting their 
rotational phase curves with low-order multipolar field models 
( [Landstreet & Mathy^[2000t [Bagnulo et al.[ 2002[ ). Although 
this techni que misses some small-scale surface magnetic field 
structures ( [Kochukhov et al.[[2004t [Kochukhov & Wad^[2010[ ), 
it is straightforward to apply, computationally inexpensive and 
therefore is widely used for obtaining information on the global 
magnetic field geometries of large stellar samples (e.g. [Auriere[ 
et al.|2007t[Hubrig et al.|2007l ). 


jects. However, this Zeeman broadening analysis is limited to 
slow rotators and provides little information about the magnetic 
field geometry. On the other hand, analysis of circular (and more 
recently linear) polarisation in spectral lines enables detection 
of much weaker magnetic fields ( [Wade et ar][200Qt [Petit et aL] 
[2011[[Marsden et al.|2014[ ), especially when polarisation signals 
can be enhanced with line-addition methods ( [Donati et al.|1997t 
[Kochukhov et al.|2010[ ). Reconstruction of the detailed surface 
magnetic field vector maps with the Zeeman Doppler imaging 
(ZDI) technique heavily relies on the interpretation of polari¬ 
sation in spectral line profiles ([Brown et al.[[ 19911 [Piskunov"^ 


[Kochukhov|2002|[Carroll et al.|2012[ ). 

For early-type magnetic stars, which typically host strong, 
globally organised, dipolar-like fossil fields, one frequently ob¬ 
serves morphologically simple, e.g. two-lobe, S-shaped, circu- 


In contrast to stable and usually topologically simple mag¬ 
netic fields of early-type stars, the dynamo-degenerated surface 
fields of active late-type stars are weak, complex and rapidly 
evolving. Their circular polarisation profile shapes are corre¬ 
spondingly more complex, often exhibiting many lobes. Such 
Stokes V profiles cannot be meaningfully characterised with in¬ 
tegral magnetic observables because their low-order moments 
tend to be close to zero even when strong polarisation signatures 
are evident in the data (e.g. [Kochukhov et al.[[2013[ ). In other 
words, the surface magnetic field distribution is so complex that 
vector averaging over the visible stellar disk leads to a substan¬ 
tial cancellation of the Stokes V signal corresponding to different 
field polarities. In this situation one usually resorts to an elabo¬ 
rate ZDI modelling of the Stokes V profiles themselves, which 
imposes stringent requirements on the signal-to-noise ratio and 
rotational phase coverage of the observational data. Being an 
intrinsically ill-posed inversion problem, reconstruction of the 
stellar magnetic field topologies with ZDI suffers from a num- 
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ber of uniqueness and reliability issues ([Donati & Brown|1997 
[Kochukhov & Piskunov||20Q2[ [Rosen & Kochukhov|2012| ), not 

the least related to the fact that only circular polarisation rather 
than the full Stokes vector spectra are typically available for ac¬ 
tive late-type stars. In this respect, ZDI is often perceived to be 
less reliable in comparison to the usual mapping of temperature 
spots or chemical inhomogeneities using intensity spectra. Due 
to a complex response of the circular polarisation spectra to the 
strength and orientation of the local magnetic field one cannot 
establish a straightforward connection between the polarisation 
profile variability pattern and major surface magnetic features in 
the same way as, for example, one can reco gnise individual sta r 
spots in the dynamic intensity spectra (e.g. [Barnes et al.|2000| ). 
It is therefore of great interest to devise simple, yet informative, 
alternative methods of the analysis of complex Stokes V profiles 
that can be used to validate ZDI results or applied when the ro¬ 
tational phase coverage is insufficient for a magnetic inversion. 

Several recent studies proposed new methods of extracting 
information from the circular polarisation signatures. [Carroll 


& Strassmeier[ ( [2014[ ) considered diagnostic potential of the net 


X +l Vl-Jc2 

dx I Ic(x,y)dy 

1 J-Vl-x2 


rotation is determined by the projected rotational velocity Ve sin i 




doVe sin i 


( 4 ) 


Under the weak-field approximation (e.g. [Landi[ 
[Degl’Innocenti & Landolfi| [2004[ ) the local intensity profile 
is unaffected by the magnetic field while the Stokes V profile is 
determined by the product of the first derivative of Stokes I and 
the line-of-sight component B^ix^y) of the local magnetic field 
vector 


V = -CzgB^(x,y)^. 

oA 


( 5 ) 


In this expression g denotes the effective Lande factor of a spec¬ 
tral line and 


Cz 


= 4.6686 X 

47rmeC^ 


( 6 ) 


absolute Stokes V signal, showing that under certain assump¬ 
tions it allows to characterise the apparent absolute longitudinal 
magnetic field. On the other hand, [Gayley & Qwockij ( [2015[ ) 
suggested to use an anti-derivative of the Stokes V profile for a 
more intuitive interpretation of the circular polarisation data. In 
this paper we further develop the latter idea. The main goal of 
this work is to present a comprehensive theoretical formulation 
of the new polarisation observable, demonstrate its connection to 
the underlying stellar surface magnetic field structure and to ap¬ 
ply the new magnetic field diagnostic method both to simulated 
circular polarisation data and to real spectropolarimetric obser¬ 
vations of stars with topologically complex magnetic fields. 


2. Cumulative Stokes V profiles 

2.1. Theoretical basis 

To describe the disk-integrated Stokes parameter profiles of a 
rotating star we consider a coordinate system with the z-axis di¬ 
rected towards the observer and the y-axis located in the plane 
formed by the line-of-sight and the stellar rotational axis. The 
star is assumed to have a unit radius and to rotate counterclock¬ 
wise as seen from the visible rotational pole. Then, the Doppler 
shift across the stellar disk is a function of the x-coordinate 
alone. With these conventions the disk-integrated continuum in¬ 
tensity is 


for the field measured in G and wavelength in A. 

Using this approximation of the local Stokes V profile one 
can express the disk-integrated circular polarisation spectrum as 


- -Czg 


d 


f>£ 






Bz{x,y) 


X — /[x,y; A-Aq- AARx]dy. 
oA 


( 7 ) 


Recalling that the observed Stokes parameter profiles are nor¬ 
malised by the disk-integrated continuum intensity, we define 
the normalised disk-integrated spectral line intensity 






X +l 

dx I I[x, y;A-Ao- AARX^dy 

1 J-y 


-VTT2 


X +l Vl-Jc2 

dx I lc{x,y)dy 

1 d- Vl-x^ 


( 8 ) 


and circular polarisation 


( 1 ) 


1 r+ Vl-jc2 

lc{x,y)dy 

-VTT? 

1 /^+ VT-T? Q 

Bz{x,y)—I{x,Y, A-Ao- AAjfXjdy. 
dA 


£'-X 


( 9 ) 


and the disk-integrated Stokes profiles of a spectral line with the 
central wavelength To are represented by the integrals 


The latter equation can be equivalently written as 


1 /^+ Vl-Jc2 



1 

Vi-x^ 

•Rv 


‘'"j 

/[x, y;A- Aq - AArx; B(x, y)] Jy 

-VT^ 

(2) 

and 


P+i 

p+ y/i-x^ 

or 

rv = 

1 dx 

I V[x, y;A- Aq - AArx; B(x, y)] Jy, 

(3) 


Ul o 

l-VT^ 

Rv 


mry-i 




Bz(x,y) 


X (Ic(x, y) - I[x, y;A-Ao- dy 


( 10 ) 


d r r' 

lly = — lCzg dx I Bz{x, 

dA y J_i J-vr:? 


where /c(x,y) represents the local continuum intensity and the 
local line intensity and polarisation are given by /(v, y, T, B) and 
y(v, y, T, B). The amplitude of the Doppler shift due to the stellar 
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O. Kochukhov: Cumulative Stokes V profiles 
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Fig. 1. Model Stokes V (upper curves) and cumulative Stokes V (lower curves) profiles for a star with several large circular magnetic spots of 
different polarities. The line profiles and the spherical maps of the radial magnetic field component are shown for three different rotational phases. 
The phases are indicated above each line profile panel. 


where ic denotes the local normalised continuum intensity 


ic(x,y) 


Ic(x,y) 

Tc 


( 12 ) 


and r/ corresponds to the local normalised residual Stokes / pro¬ 
file 


ri[x,y\A - Aq - I^Arx] 


Ic(x, y) - I[x, y;A-Ao- AArx] 
Ic(x,y) 


( 13 ) 


We define the cumulative Stokes V (CSV) profile as 


QvW = 


1 

Wi 



nv(A')dA\ 


( 14 ) 


where 


Wi = 



(1 - ni)dA 


( 15 ) 


is the equivalent width of the disk-integrated Stokes I profile and 
the integration limits Tmin, ^max cover the full extent of a spectral 
line. Substituting Eq. ^TT) into Eq. ( p^ yields 

CyP ^+VT^ 

Wi 

X ic{x,y)ri[x,y\ A-Aq- AARx]dy. (16) 


Dividing this quantity by Czg provides the longitudinal mag¬ 
netic field density 


BfiA) ^ 


@v(^) 

Czg 




Bz{x,y) 


X ic{x,y)ri{x,y\ A - Aq - AARx\dy. 


( 17 ) 


This quantity represents a velocity-resolved measure of the line- 
of-sight magnetic field component, weighted by the projected 


surface area and by the equivalent width-normalised local Stokes 
I profile. Bz{A) has the units of magnetic field divided by wave¬ 
length. Its integration over the full line profile gives the nor¬ 
malised first moment of Stokes V 



B,{A)dA = - 


1 

CzgWi 



(A-AoWA^{B,}, 


( 18 ) 


commonly known as the mean longitudinal magnetic field. 

The QvW observable discussed above is identical to the 
Stokes V anti-derivative concept introduced by [Gayley & 


Owockij ( |2015| ) with the exception that we formulate this quan¬ 
tity with respect to the normalised Stokes V parameter "Ry and 
divide the resulting profiles by Wi whereas the earlier study pro¬ 


posed to divide hy 1- Ri. The quantity Qy(A) discussed by Gay- 
|ley & Qwocki| ( |20f5] ) is easily recovered from our CSV profiles 
using the transformation 


V- Ki 

The corresponding velocity-resolved mean longitudinal mag¬ 
netic field B^(A), which can be obtained by dividing ^y(A) by 
Czg, has the units of magnetic field strength. 

The longitudinal magnetic field density B^iA) or the velocity- 
resolved longitudinal field B^ (A) derived from the CSV profiles 
represent morphologically simpler observables compared to the 
Stokes V profiles themselves. They characterise the line-of-sight 
magnetic field averaged along the stripes of constant Doppler 
shift, thus providing a more intuitive representation of the infor¬ 
mation content of the Stokes V spectra. 

Although the observables B^iA) and B^ (A) are closely related, 
different normalisation choices lead to a somewhat different be¬ 
haviour. The velocity-resolved mean longitudinal field closely 
tracks the actual local line-of-sight magnetic field component. 
The longitudinal field density includes an additional weighting 
by the projected surface area and therefore gradually diminishes 
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to zero towards the line profile edges even if the underlying 
distribution is uniform. The (T) quantity compensates chang¬ 
ing projected area by the \-Ki normalisation but instead suffers 
from major noise artefacts at profile edges due to division of two 
small numbers. Compared to B^iX), this quantity may also be 
adversely affected by the distortions of Stokes I profile shape 
caused by stellar surface inhomogeneities. 

Equation ^TJ) can be considered as a convolution integral 
with a kernel given by the residual Stokes I profile. The larger 
is the ratio of AAr to the local profile width AT/, the less suscep¬ 
tible is Bz(A) to the cancellation of polarisation signatures cor¬ 
responding to different field polarities. Therefore, the new mag¬ 
netic observable will be most effective for rapid rotators since in 
this case AAr » AT/. 

Based on thes^e definitions, we carried out an illustrative cal¬ 
culation of the ^^(T) profiles for a magnetic field distribution 
comprising several circular spots with the radial field strength 
of 0.5-1.0 kG (see Fig. [^. The calculations assumed a Gaus¬ 
sian shape with FWHM = 8 kms“^ for the local Stokes I pro¬ 
file. The adopted central wavelength and effective Lande factor. 
To = 5630 A and g = 1.215, corresponded to the mean values of 
the LSD line mask used for V410 Tau (see Sect.|^. The star was 
assumed to rotate with Ve sin i = 50 km s“^ and to be inclined by 
i = 60° with respect to the line of sight. 

Figure shows the spherical maps of the radial magnetic 
field as well as the Stokes V and CSV profile for three different 
rotational phases. When the surface magnetic field distribution 
is dominated by one of the magnetic polarities, the Stokes V 
profile exhibits the classical S-shape pattern. The corresponding 
CSV profiles show a single bump, which can be either positive 
or negative, depending on the dominant magnetic field polar¬ 
ity. On the other hand, when several regions of different polarity 
are present on the stellar disk (e.g. middle panel in Fig. [^, the 
Stokes y profile exhibits multiple components, which usually 
cannot be interpreted in terms of the surface magnetic field dis¬ 
tribution without carrying out a detailed modelling. In contrast, 
the CSV profiles for the same rotational phase are simpler and 
readily show spots of different polarity at the longitudes corre¬ 
sponding to the Doppler shift within a spectral line profile. 


2.2. Deconvolved longitudinal magnetic field density 

A simplified form of Eq. can be obtained by assuming that 
the local residual intensity profile is represented by a Gaussian 
function which does not vary across the stellar surface. This ap¬ 
proximation, together with the weak-field assumption, is com¬ 
monly made in the context of modelling stellar circular polari¬ 
sation spectra (e.g. [Petit et al.|2004l[Petit & Wade|2012| ). Then, 
omitting To, the normalised residual profile becomes 


ri(A, x) 
Wi 


1 




exp 


(A - AA/,xf 
2cr^ 


( 20 ) 


Since the kernel described by this equation is known once cr and 
Ve sin i are specified, it is possible to define deconvolved longitu¬ 
dinal magnetic field density 


X + Vl-x^ 

Bz{x,y)ic{x,y)dy, 

VTT2 


( 21 ) 


in which the effect of averaging over the line profile width is 
taken out. 
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Fig. 2. Theoretical dynamic Stokes V (left panel) and cumulative 
Stokes V (right panel) profiles corresponding to the radial magnetic field 
maps shown in Fig.[^ 


In practice, ^*(T) can be obtained from the observed BfiA) 
with the help of a suitable deconvolution algorithm. For the in¬ 
terpretation of B^ (T) it can also be assumed that the centre-to- 
limb variation of the continuum intensity is described by a linear 
limb-darkening law 


icix,y) 


3(1 - s + sji) 
7r(3 - d) 


( 22 ) 


where £ is a limb-darkening coefficient and ju = - x'^ -y'^. 


2.3. Dynamic CSV profiles 

A two-dimensional plot of the Steaks I profile variability pattern 
as a function of wavelength (or velocity) and rotational phase 
is commonly used to assess distribution of the stellar surface 
inhomogeneities. Although similar plots of the Stokes V pro- 
file va riation are also occasionally published ( [Donati et al.|1999 


[2006[ ), they could not be directly employed for obtaining infor¬ 
mation about the stellar surface magnetic field distributions. The 
dynamic cumulative Stokes V profiles, on the other hand, allow 
one to trace individual magnetic spots and characterise the stellar 
magnetic field geometry. 

An example of the dynamic CSV profiles is shown in Fig.j^ 
This figure presents polarisation spectra and corresponding B^ 
profiles for the model magnetic field distribution discussed in 
Sect. [2.1[ The longitudinal positions and polarities of the four 
magnetic spots can be directly read out from the dynamic CSV 


























O. Kochukhov: Cumulative Stokes V profiles 
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Fig. 3. Effect of a random noise with cr = 2 X 10 on the Stokes V 
{upper panel) and CSV profiles with Wi {middle panel) and \-Ki {bot¬ 
tom panel) normalisations. The thick double line shows the noise-free 
calculations corresponding to the rota tional phase 0.125 of the model 
Stokes profiles discussed in Sect. |2.1| The thin lines represent profiles 
for different noise realisations. The and profiles are obtained with 
the forward integration according to Eq. ( p^ . 


plot. The latitudes of the spots can be deduced from the velocity 
span of their signatures in the profile. 

2.4. Deriving CSV profiles from noisy data 

The cumulative Stokes V profiles exhibit non-trivial noise prop¬ 
erties. On the one hand, since the CSV spectra are computed by 
integrating the Stokes V signatures, one may expect some can¬ 
cellation of the random noise. On the other hand, the resulting 
noise in the CSV profiles themselves is highly correlated, lead¬ 
ing to a qualitatively different behaviour compared to the ini¬ 
tial Stokes V profiles. These special noise properties have to be 
considered when applying CSV diagnostic to real observational 
data. 

We investigated the noise properties of the CSV profiles with 
several sets of Monte-Carlo simulations. First, we considered the 




Fig. 4. Upper panel: cumulative distribution of xl values for a set 
of independent, normally-distributed random variables (solid line) and 
foxxl corresponding to Qv{^) given by Eq. (24i (dashed line). Lower 
panel: fraction of detections of magnetic fielo^gnatures in the Stokes 
V (solid line) and CSV profiles (dashed line) as a function of random 
noise added to the model Stokes profiles. 


model Stokes V spectra corresponding to the four-spot magnetic 
field distribution discussed in Sect. |2.1| The mean peak-to-peak 
amplitude of these circular polarisation profiles is 2x10“^. The 
profiles were sampled with a step of 1 km s“\ yielding about 100 
spectral points across the line. A random, normally-distributed 
noise with cr = 2 x 10“^ was added to these profiles. Typical 
CSV signatures resulting from the forward integration according 
to Eq. ( p^ are presented in Fig. It is evident that, while the 
CSV spectra appear to have less random noise compared to the 
initial Stokes V profiles, they suffer from a systematic, ramping 
deviation from the zero line, which increases as the integration 
proceeds from blue to red. The bottom panel of Fig. [^illustrates 
the CSV observable normalised by the residual profile inten¬ 
sity 1 - ^/ according to Eq. ( [T9| ). The devastating impact of the 
noise at profile edges is apparent. 

It is possible to define the CSV profiles using an alternative, 
backward, red-to-blue integration scheme 

= M')dA'. (23) 

Ja 

This formula gives identical results to Eq. ( [T^ in the absence of 
noise. When the noise is present, it will lead to a systematic dis¬ 
crepancy similar to the one illustrated in Fig. but increasing 
towards the blue side of the line. The overall systematic devi¬ 
ation from zero can be minimised by a weighted mean of the 
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forward and backward integrations 


_ 1 r W Wmax 

I Ky{A')dA' —w I *Rvi^')dA' 

where the weights w- are given by 

+ d-iYiax ~ A _ A — /Ijjiin 

w =- and w = --. 

-^max “ -^max -^max “ -^max 


(24) 


(25) 


The correction expressed by these equations is equivalent to sub¬ 
tracting from the original CSV profile Qv(^) a straight line going 
through 0 at /Imin and ^y(/ljjiax) at /ljjiax* 

Application of Eqs. ( [24}j25] ) still results in a substantially 
higher reduced than found for the original Stokes V pro¬ 
files corrupted by noise. Using 10^ realisations of normally- 
distributed random noise in the absence of any background sig¬ 
nal, we have built a cumulative distribution function for the x^ 
values deduced by performing integration according to Eq. p4| ) 
with associated error propagation. As shown in the upper panel 
of Eig.|^ integration leads to substantially largervalues com¬ 
pared to the initial Stokes V profiles. Eor example, the 10“^ 
probability threshold corresponds toXy ^ 6 for the CSV profiles 
compared toXy ~ 1.4 for the initial profiles. 

Based on this empirical CDE ofXy of CSV profiles, we esti¬ 
mated the fraction of such profiles that will be detected above the 
noise level with the false alarm probability of 10“^ or smaller for 
the four-spot magnetic field geometry sampled at 100 equidis¬ 
tant rotational phases. These MC simulations were carried out 
for noise amplitudes in the range from 2 x 10“"^ to 2 x 10“^, i.e. 
10-100% of the mean peak-to-peak Stokes V profile amplitude. 
Simulations were repeated 10"^ times for each value of the noise 
amplitude. 

As can be seen from the lower panel of Eig. the detection 
of magnetic field signatures using CSV profiles is less successful 
than with the original Stokes V spectra. Therefore, the CSV 
diagnostic technique does not represent a useful alternative, in 
terms of the mere field detection, to the widely used -based 
Stokes V profile detection method. 


3. Observational data 


Eor the purpose of testing the CSV diagnostic method on real 
data we analysed high-resolution spectropolarimetric observa¬ 
tions available from the Polarbase archive ( [Petit et ar]|2014| ). 
This database provides reduced, calibrated, one-dimensional 
Stokes parameter spectra obtained with ESPaDOnS at CEHT 
and Narval at TBL of the Pic du Midi observatory. Both in¬ 
struments are fibre-fed, thermally stabilised, echelle spectropo- 


larimeters covering the wavelength range from 3694 to 10483 A 
in a single exposure at the resolving power of A/AA = 65 000. 

As an example of the Stokes V profiles of a cool active star 
with a complex magnetic field we considered spectropolarimet¬ 
ric observations of the weak-line T Tauri star V410 Tau. This 
data set, dis cussed in the papers by jSkelly et al.| ( |MT0] ) and jRicej 
et al.| ( |20lT] ), comprises 45 Stokes V observations obtained in the 


span of about two weeks in January 2009. The signal-to-noise 
ratio of these spectra, S/N 150, is insufficient for the detec¬ 
tion and analysis of circular polarisation signatures in individual 
spectral lines. Therefore, we applied the least-squares decon¬ 
volution approach ( jKochukhov et al.||2010| ) to each spectrum, 
combining about 4400 spectral lines into a single high S/N ratio 
profile. The data were phased with the ephemeris ofjStelzer et al.j 
( [20031 ). 


We also analysed the Stokes V observations of the He-strong 
Bp star HD 37776. This object is known to possess a complex 
magnetic field, strongly deviating from an oblique dipolar ge¬ 
ometry ( jKochukhov et al.|[2011| ). There are 27 Stokes V spec¬ 
tra of HD 37776, with a typical S/N ratio of 450, available in 
Polarbase. These data were acquired over the period of 2006- 
2012. We used the variable-period ephemeris of jMikulasek et al.j 
( |2008| ) to compute rotational phases. Since the magnetic field 
of HD 37776 is very strong, circular polarisation signatures are 
readily observable in individual spectral lines. In this paper we 
study the He i 6678 A line. Its circular polarisation shows a vari¬ 


ability pattern very similar to the He i 5876 A line analysed by 
jKochukhov et al.| ( '2011| ) using a lower quality data. 


4. Examples of the CSV diagnostic 

4.1. Weak-line T Tau star V410 Tau 

The young rapidly rotating star V410 Tau (HD 283518) exhibits 
ample signs of the surface magnetic activity and was targeted by 
a number of temperature DI studies ( jHatzesjl 995} [Rice & Strass- 
|meier|1996[|Rice et al.|201lj ). These analyses revealed a persis¬ 
tent large polar spot and evolving temperature inhomogeneities 
at lower latitudes. Magnetic field was detected on V410 Tau by 
Donati etar] ( |1997| ). jSkelly et ah ( |201Q| ) reconstructed the sur¬ 


face magnetic field topology of this star with ZDI, finding a sig¬ 
nificant toroidal magnetic component and a complex radial field 
distribution. On the other hand, [Carroll et al.j ( [2012[ ) reported a 
relatively simple, predominantly poloidal magnetic field struc¬ 
ture consisting of two large polar magnetic spots with an op¬ 
posite polarity. The Stokes V spectropolarimetric observations 
interpreted by jSkelly et al.[ ( [MT0| ) and jCarroll et al.[ ( [2012[ ) were 
obtained with the twin instruments (ESPaDOnS and Narval) and 
at practically the same time, meaning that the discrepant mag¬ 
netic inversion results cannot be ascribed to the intrinsic varia¬ 
tion of the surface magnetic field of V410 Tau. 

Here we use the combined ESPaDOnS and Narval spec¬ 
tropolarimetric data set to reconstruct independent maps of the 
magnetic field topology and brightness distribution. The inver¬ 
sion methodology applied to V410 Tau is described in detail by 
jKochukhov et al.[ ( [2014[ ). Briefly, the magnetic field of the star 
is represented in terms of a superposition of the poloidal and 
toroidal harmonic terms with the angular degree up to ^ = 15. A 
penalty function prohibits unnecessary contribution of the high- 
order harmonic modes. A separate regularisation procedure ap¬ 
plied to the brightness map minimises any deviation from the 
reference (photospheric) brightness value. We employed the an- 
alytical Unno-Rachkovsky Stokes parameter profiles (e.g. jUandT 
jDeglTnnocenti & Landolfij[2004j ) to approximate the local in¬ 
tensity and circular polarisation spectra. The brightness and 
magnetic field maps were recovered self-consistently from the 
LSD Stokes I and V profiles, adopting i = 60° and Vg sin i = 
74 kms“^ 

The LSD profile fits and the resulting brightness and mag¬ 
netic field maps are presented in Eig.[^ The surface distributions 
which we obtain are generally compatible with the maps pub¬ 
lished by jSkelly et al. ( j2010j ). Similar to these authors we find 
a dominant dark polar spot and numerous small-scale brightness 
inhomogeneities at lower latitudes. The poloidal and toroidal 
field components contribute 63% and 37% respectively to the 
total magnetic field energy. The field strength reaches ~ 1 kG 
locally on the stellar surface. The azimuthal field exhibits two 
large unipolar regions while the radial field shows a more com¬ 
plex structure. The darkest areas in the brightness map do not co- 
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Fig. 5. Results of the ZDI analysis of V410 Tau. The four panels on the left side of the figure compare the observed (histogram) and theoretical 
(solid line) LSD Stokes I and V profiles. In these plots the spectra corresponding to different rotational phases are shifted vertically. The phase 
is indicated to the right of each profile. Reconstructed maps of the radial, meridional and azimuthal magnetic field components as well as the 
brightness distribution are presented in the right column. The star is shown using the flattened polar projection between latitudes -60° and +90°. 
The thick circle corresponds to the rotational equator. The colour bars give the field strength in kG and the brightness relative to the photospheric 
value. 


incide with any particular magnetic field features. We have also 
verified that a very similar magnetic field geometry is obtained if 
one performs a direct ZDI reconstruction of each magnetic field 
component without using the spherical harmonic formalism. 

Based on the magnetic inversion results, we calculated the 
CSV profiles using both the observed and the ZDI model Stokes 
I and V spectra of V4I0 Tau. Figure compares profiles ob¬ 
tained by applying Eq. ([g, Eq. ( |^ , and the weighted mean of 
the forward and backward integrations given by Eq. ( [24| ). It is 
clear that, in the first and second case (Figs. |^b), a significant 
systematic deviation between the observed and theoretical CSV 
profiles appears due to the noise accumulation. On the other 
hand, systematic effects are largely brought under control in the 
third set of CSV profiles (Fig.[^). 


Figure illustrates the impact of employing the I-"Rj nor¬ 
malisation in place of the Wj normalisation used elsewhere. The 
profiles shown in this panel directly provide the average line- 
of-sight magnetic field component at different longitudes on the 
stellar surface. However, these profiles, being morphologically 
similar to the spectra, are visibly distorted by the noise arte¬ 
facts at the line edges. 

We tested the accuracy of the CSV diagnostic with the model 
LSD profiles corresponding to the magnetic field map shown in 
Fig. For clarity we disregarded the non-uniform brightness 
distribution and considered the line profiles for 10 equidistant 
rotational phases. The CSV profiles were computed from the 
simulated observations using Eq. ( p^ and by performing a nu¬ 
merical integration of the line-of-sight magnetic field component 
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Fig. 6. Cumulative Stokes V profiles of V410 Tau computed from the observed (histogram) and theoretical (thick solid line) LSD circular 
polarisation spectra of V410 Tau shown in Fig. The CSV profiles are shifted vertically according to their rotational phases similar to the 
Stokes V profiles in Fig.[^ The first three panels show the CSV spectra computed with a) forward integration, b) backward integration, and c) 
weighted mean of the backward and forward integration. The last panel d) shows the result of converting the CSV profiles in panel c) to the 1 - “R/ 
normalisation according to Eq. ( p^ . 


according to the right-hand side of Eq. ^Tf) . In these calculations 
we used the same local non-magnetic rj profile as was adopted 
for the ZDI inversions. As expected, the two sets of CSV profiles 
agree perfectly (Fig. [7^). 

The CSV method is based on the weak-field approximation, 
which is usually considered to be valid only for the magnetic 
field strengths below ~ 1 kG. It is useful to assess the usefulness 
of the cumulative Stokes V profiles for much stronger magnetic 
fields. To this end, we scaled all the vector components of the 
ZDI map of V410 Tau by a factor of 10. This increased the mean 
field strength from 570 G to 5.7 kG. We then repeated the direct 
and geometrical evaluation of the CSV profiles. The resulting 
spectra are compared in Fig. [^. The two sets of profiles 
show only marginal discrepancies despite that the characteristic 
magnetic field strength now significantly exceeds 1 kG. There¬ 
fore, it appears that the CSV diagnostic is usable well beyond 
the nominal limit of the weak-field approximation. 

Using the ZDI maps of V410 Tau we calculated variation 
of the Stokes parameter profiles with a dense rotational phase 


sampling. The resulting dynamic circular polarisation and CSV 
spectra are presented in Fig.[^ Unlike Fig.j^it shows a stationary 
pattern (phases 0.0-0.3) in addition to features travelling across 
the stellar disk (e.g. phase 0.6-0.8). The travelling features 
are associated with the radial field spots, in particular the neg¬ 
ative polarity spot best visible at the rotational phase 0.75. The 
stationary features correspond to the large unipolar azimuthal 
fields regions. Such behaviour of the dynamic CSV profiles can 
be considered as an evidence for the presence of a significant 
toroidal magnetic field on the stellar surface. 

4.2. He-strong Bp star HD 37776 

The B-type magnetic chemically peculiar star HD 37776 (V901 
Ori) is known for its unusually complex, double-wave longitu¬ 
dinal magnetic field variation ( [Thompson & Landstreet|p^85] ). 
This behaviour indicates a large deviation of the stellar sur¬ 
face magnetic field topology from an oblique dipolar geome¬ 
try. In the literature this object is often considered to be an 
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O. Kochukhov: Cumulative Stokes V profiles 



Fig. 7. Theoretical Stokes /, V and CSV profiles corresponding to a) the magnetic field map of V410 Tau shown in Fig.|^and b) the same map 
scaled by a factor of 10. In each case the CSV column compares 5^(4) obtained from the Stokes I and V line profiles (solid line) with the results 
of numerical integration of the line-of-sight magnetic field component over the stellar disk (dashed line). 


example of the star with a quadrupolar magnetic field struc¬ 
ture. Indeed, |BohIender| ( |1994| ) and [Khokhlova et al.| ( |2000| ) 
derived quadrupolar magnetic geometry models for HD 37776 
based on fitting the longitudinal field curve and Stokes V pro¬ 
file variation, respectively. However, using a more extensive set 
of phase-reso lved spectropolarimetric observations, [Kochukho^ 
jet al.j ( |201I| ) showed that the circular polarisation spectra of 
HD 37776 cannot be fitted with an axisymmetric quadrupolar 
field structure. Instead, the ZDI inversion carried out in that 
study suggested a complex field distribution comprising a series 
of magnetic spots with alternating positive and negative polarity. 


The high-resolution circular polarisation spectra of 
HD 37776 analysed here enable a straightforward verifica¬ 
tion of the ZDI results of [Kochukhov et al.j ( [2011[ ) using the 
cumulative Stokes V profile diagnostic method. Figure shows 
the Stokes V and CSV spectra centred on the Hei 6678 A 
line. The dynamic CSV profiles are dominated by a series of 
travelling features. There is no evidence of stationary structures 
like in the case of V410 Tau, indicating that no strong toroidal 
magnetic field is present. Furthermore, a careful examination 
of the CSV profiles allows one to identify three spots with 
negative polarity and the same number of positive-field features. 
This corresponds very well to the six well-defined spots seen 
in the radial magnetic field map recovered with ZDI (see fig. 5 
[Kochukhov et al.[[20II[ ). On the other hand, considering the 
dynamic CSV profiles in Fig. one can immediately rule out 
the quadrupolar field hypothesis. A quadrupolar magnetic field 
geometry with a large obliquity, as required to reproduce the 
(5z) variation of HD 37776, exhibits four regions of alternating 
field polarity along the stellar rotational equator. But Fig. 
shows at least six regions. 


5. Summary and discussion 

In this paper we examined the information content of the cumu¬ 
lative circular polarisation profiles and assessed their usefulness 
for the analysis of stellar surface magnetic fields. This work rep¬ 
resents a further development of the idea put forward by [Gayle^ 
[& Qwocla| ( [2015| ). Starting with the weak-field approximation, 
we gave a definition of the new magnetic observable as a nor¬ 
malised integral of the Stokes V spectrum and showed its relation 
to the underlying surface magnetic field structure. The transfor¬ 
mation from the Stokes V to the CSV profiles greatly simplifies 
the morphology of polarisation spectra and enables a direct and 
intuitive identification of the regions of different field polarity on 
the stellar surface. The CSV observable characterises the line- 
of-sight magnetic field component weighted by the continuum 
intensity and by the Doppler-shifted local residual Stokes I pro¬ 
file. The CSV spectrum thus provides a velocity-resolved mea¬ 
sure of the stellar longitudinal magnetic field. Consequently, this 
observable is most useful for rapidly rotating stars with complex 
magnetic field topologies. 

Using theoretical circular polarisation profile calculations 
and real observational data we assessed visibility of the sur¬ 
face magnetic field features in the dynamic CSV profiles. Un¬ 
like the dynamic Stokes V spectra, which can hardly be inter¬ 
preted directly, the dynamic CSV profiles enable a qualitative 
analysis of the magnetic star spot distributions. With the help 
of Monte Carlo simulations we studied the effect of a random 
observational noise on the CSV profiles. The noise signatures 
in the CSV spectra are highly correlated, leading to a system¬ 
atic ramping offset from a noise-free version of the line profiles. 
We showed that this problem can be minimised using a weighted 
mean of the CSV profiles computed with the forward and back¬ 
ward integration of the original Stokes V data. Nevertheless, the 
presence of correlated noise in the CSV profiles generally makes 
them inferior for the purpose of field detection compared to the 
original Stokes V spectra. 
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Fig. 8. Theoretical dynamic Stokes V (left panel) and cumula- Fig. 9. Observed dynamic Stokes V (left panel) and cumulative Stokes 

tive Stokes V (right panel) profiles corresponding to the ZDI map of V (right panel) profiles of the He i 6678 A line in Bp star HD 37776. 

V410 Tau shown in Fig.[^ 


With the goal of testing the CSV magnetic field diagnos¬ 
tic methodology we studied the archival high-resolution circu¬ 
lar polarisation observations of the cool active star V410 Tau 
and the hot magnetic star HD 37776. Both objects are rapid ro¬ 
tators with complex surface magnetic field topologies. For the 
weak-line T Tauri star V410 Tau we obtained brightness and 
magnetic field distributions with the help of the ZDI modelling 
of the least-squares deconvolved Stokes I and V profiles. The 
resulting magnetic map bears a close resemblance to the mag¬ 
netic field topology obtained by |Skelly et al.| ( |2010| ) from the 
same data. However, our i nversion results cont radict the mag¬ 
netic field map published by [Carroll et al.| ( [20T^ whose data set 
was also included in our analysis. Our ZDI modelling accounted 
for the effect of an inhomogeneous brightness distribution on the 
Stokes V profiles. We have also verified that very similar mag¬ 
netic field maps are obtained using the direct, pixel-based ZDI 
and the magnetic inversions relying on the spherical harmonic 
parameterisation of the stellar magnetic field. Therefore, these 
two aspects are unlikely to be responsible for the discrepancies 
of the ZDI map s presented here and in jSkehy et al.| ( |2010| ) on the 
one hand and in jCarroll et al.| ( [MT^ on the other hand. It is more 
likely that this disagreement is rooted in the differences between 
the LSD (used here and by jSkehy et al.|201()| ) and PC A (used by 
[Carroll et al.|2Q12j ) mean polarisation profiles. A comprehensive 
comparative study of the LSD and PCA line-addition techniques 
is required to address this problem. 
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The dynamic CSV profiles of V410 Tau indicate the pres¬ 
ence of a strong azimuthal magnetic field component on the stel¬ 
lar surface. This agrees with our ZDI results and suggests a 
non-negligible contribution of the toroidal magnetic field. We 
also used the ZDI map of V410 Tau to verify interpretation of 
the integral of Stokes V profiles in terms of the disk-averaged, 
weighted line-of-sight magnetic field component. These calcu¬ 
lations were repeated for a scaled-up version of the stellar mag¬ 
netic field geometry, demonstrating that the usefulness of the 
CSV diagnostic extends far beyond the nominal limit of the weak 
field approximation. 

For the magnetic Bp star HD 37776 we analysed the CSV 
dynamic spectrum of the Hei 6678 A line. The variation of its 
CSV profiles is clearly inconsistent with the quadrupolar mag¬ 
netic field topology frequently suggested for this star. On the 
other hand, multiple magnetic spots which can be identified in 
the dynamic CSV spectra of HD 37776 are consistent with the 
surface magnetic featur es in the ZDI map of this star derived by 
[Kochukhov et al.[ ( [201 1[ ). The CSV analysis thereby confirms the 
complex, non-quadrupolar nature of the magnetic field topology 
of HD 37776. 
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